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Abstract 



The large-Af behavior of the quartic-cubic generalized two dimensional 
Yang Mills V{N) on the sphere is investigated, for small cubic couplings. 
It is shown that single transition at the critical area which is present for the 

'q\ '• quartic model, is split into two transitions, both of them are third order. 

(~| I the phase diagram of the system for small cubic couplings is obtained. 



. . 1 Introduction 

^— N ' During recent years, there has been extensive studies of the two-dimensional 

^ ! Yang-Mills theory (YM2) and generalized Yang-Mills theories (gYM2's) [1-18]. 

CO ' A review on this topic is [19]. These are important integrable models which are 

r — , expected to shed light on some basic features of pure QCD in four dimensions. 

f^ ' Besides, there are certain relations between these theories and string theories. 

00 . To be more specific, there relations are between the large-gauge-group limit of 

^^ ' YM2 and string theory. For example, it is shown in [3], [6], and [7], that a gauge 

theory based on SU(iV) is split at large N into two copies of a chiral theory, 

encapsulating the geometry of the string maps. The chiral theory associated to 

the Yang-Mills theory on a two-manifold A^ is a summation over maps from 

C^ ' the two-dimensional world sheet (of arbitrary genus) to the manifold A4. This 

leads to a 1/A^ expansion for the partition function and observables, which is 
convergent for all of the values of area x coupling constant on the target space 
A4, if the genus is one or greater. 

The partition function and the expectation values of the Wilson loops of 
YM2 were obtained for theories on a lattice [1,20] and continuum [4,11,12]. 
The partition function and the expectation values of the Wilson loops of gYM2's 
were calculated in [9,10]. All of these results are in terms of summations over the 
irreducible representations of the corresponding gauge group. When the group 
is large, these summations are dominated by some specific representations (the 
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so called classical representations) in some cases, and one can obtain closed-form 
expressions for that representation and the observables of the theory. 

The large- TV limit of U(A'^) YM2 on a sphere was first studied in [21]. There 
the classical representation was calculated and it was shown that the free energy 
of theory has a logarithmic behavior with respect to the area of the sphere. In 
[13], it was shown that that behavior is correct provided the area of the sphere is 
smaller than some critical area, and it was further shown that at that area a third 
order phase transition occurs. This transition is similar to the well-known Gross- 
Witten-Wadia phase transition for the lattice two dimensional multicolor gauge 
theory [22,23]. The phase structure of the large- iV YM2, generalized YM2's, 
and nonlocal YM2's on a sphere were further discussed in [16, 18,24,25]. The 
large- A'' limit of the partition function of YM2 on orientable compact surfaces 
with boundaries was discussed in [26] , and the large- iV behavior of Wilson loops 
of YM2 and gYM2 on sphere were investigated in [27]. The critical behaviors 
of these quantities have been also studied. 

A gYM2 with the gauge group V{N) on a sphere is characterized by a func- 
tion G, as introduced in [18], for example. Most of these investigations on 
gYM2's have been on theories with even G. In fact, as it was pointed out 
in [18], if G is a polynomial then its degree should be even in order that the 
partition function be convergent. So if one takes G a monomial, it should be a 
monomial of even degree. So to study theories based on non-even G's, the corre- 
sponding G's should contain at least two terms. A quadratic-linear combination 
gives nothing new, as the linear term can be absorbed in a shifted variable. So 
the simplest nontrivial non-even G must contain a quartic term. 

In this paper the behavior of a gYM2 based on the gauge group V{N) on a 
sphere is studied, for which the function G is a binomial containing a quartic 
term as well as a cubic one. The study is performed for small values of the 
cubic coupling. In section 2 the behavior of this model for areas smaller than 
the critical area (the weak phase) is investigated and compared with the behavior 
of a theory without the cubic term. In section 3 the transitions properties at 
the first transition area are investigated. In section 4 it is seen that there exists 
another transition, and the relation between this transition to the first one is 
investigated. Section 5 is devoted to the concluding remarks. 

2 The weak phase 

To fix notation, let us first quickly review the expression for the partition func- 
tion of a \J{N)-gYM2 on the sphere, in the large N limit [18]. 
The partition function of a gYM2 on the sphere of area A is 

Z = Y^dle-^^(^\ (1) 

r 

where r's label the irreducible representations of the gauge group, dr is the 
dimension of the representation r, and A is a combination of the Casimirs of 



the group. For the gauge group IJ{N), one can take A as: 

Mr) = tj^C,ir), (2) 

k=l 

in which Ck is the fc'th Casimir of the group, and a^'s are arbitrary constants. 
The representations of the group V{N) are parametrized by N integers ni > 
n2 > ... > n^. In terms of this parametrization, one has 



dr= n (i 



l<i<j<N ^ ■'' 

N 

Ck ^Y.i(n^ + N - i)'' - {N - i)% (3) 

i=l 

In order that the partition function (1) be convergent, it is necessary that p in 
Eq. (2) be even and ap be positive. 

For the case of large N, the summation can be rewritten in the form of a 
path integral over continuous parameters. These parameters arc introduced as 

< a; := — < 1, 

n{x) := -, 

(j){x) := - n{x) - 1 + X. (4) 

The partition function (1) then becomes 

Z^ f H d<f,{x) exp[5(^)], (5) 

•' 0<2;<1 

where 

S{c^) -.^ N^ La 1^ dxG[0(x)]+^ dxj^ dy log|0(x)-0(y)||, (6) 
apart from an unimportant constant, and 

G(z):=f](-l)'=a,z^ (7) 



fc=i 



Introducing the density 



one would have 



M*(.^)l - ally. (8) 



dz p(z) = 1, (9) 



(the so called normalization condition for the density), where b and a are the 
lower and upper limits of 4>{x), respectively. The condition that n^'s be nonin- 
creasing with respect to i imposes the following condition on the density. 

< p{z) < 1. (10) 

The path integral in the right-hand side of (5) is dominated by the classical 
representation, which maximizes the function S, or minimizes the free energy 
F defined through 

F:=-^\nZ. (11) 

Denoting the density corresponding to this representation by pd, it is seen that 
the derivative of the free energy with respect to the area is 

F'{A) ^ f dx G[0ci(x)], 

JO 



dz pci{z)G{z), 

f dz 
= f :^Hdiz)Giz), (12) 

Jc^ 2 7ri 

where the function iJ of a complex variable is defined through 

H{z):= TdC^, (13) 

Jb z — (, 

and Coo is a counterclockwise contour outside of which H is analytic. It is seen 
from (9) that H{z) behaves like z~^ for large z. 

The classical density pd is the weak density p^ satisfying 

g{z) = P r dC ^^, b^<z< aw, (14) 

where P means principal value and 

giz):=^G'{z), (15) 

provided the conditions (10) are not violated by /?„. This is the weak phase. 
One obtains 



Hw{z) = g{z) - \J(z- Ow) (z - 6w) 

X V (2»~l)n(29-l)!! „ („+^+,+r)(o) (^g) 

A^ 2"+9n!(7! n + g + TO+1 ! ^^ ^ \ )^ \ t 

where g^"-* is the n'th derivative of g, and 



, . y/{aw - z){z~ 5w) 
Pw(z) = 



TT 



E (2^-l)!!(2g-l)!! alblz-g^-^-^^+'Ho). (17) 

Z^ 2^+1 n\ q\ (n + q + m + ly. ^ ^ -^ \ ) K / 

m,n,q—0 a \ a / 



The condition that H{z) behaves Hke z ^ for large z, is equivalent to the fol- 
lowing equations. 

- (2n l)n(2,-l)M 

These equations are used to obtain a^ and fe^. 
Now consider a function G like 



G{z)^z^ + Xz\ 




(20) 


The conditions (18) and (19) become 






/ 3 \ / 3 \ 

T^ (ScTw + ^AJ +Cr^ (2(Tw + -AJ 


= 0, 


(21) 


^Tw + Tw ( 3cr^ + -crwAj 


1 


(22) 


where 

a + & 
(T := , 







2 
a — b 



Using (16), one obtains 



' w 



(3 A) 



(23) 



^w(^) -^{4z3 + (3z2)A_^(z_a„)2_^2 

X [44 + 2r2+4awZ + 4z2 + 3((T„ + z)A]}. (24) 

The solution to (21) and (22) for small values of A is 

1 (3A)V2 3 
aw= -^A ^^— A +0(A ), 



_ + A^ + 0(A4). (25) 



Using these, (12), and (24), one obtains a perturbative expression for the deriva- 
tive of the free energy with respect to the area: 

To investigate the conditions (10) for p^, one obtains p^ and its minima and 
maxima. From (17), one has 

Pw{z) = -- x/t^ -{z- crw)2 [4cr2 + 2t2 + 4crw z + 4z2 + 3 (cTw + ^) A]. (27) 

Z TT 



Using this, one obtains three points where the derivative of pw vanishes: 



and the corresponding values for the density p^: 

Pw{zo) 



[3 A) 



1/2 



96 



^<"''" <"'•'" A= + 0(A'), 



A3 + 0(a4), (28) 



Stt 



a287r 



p„(.±) = i%il^T-^A3 + 0(A'^). 



Stt 



16 TT 



(29) 



The absolute maximum of pw exceeds one at A = Ac. This gives the critical 
area -/in ■ 



Ac=Ae(0) 



1- 



^^^°)|A|3 + 0(|An 



47r 



where 



^c(O) 



27 TT^ 

256 



(30) 
(31) 



It is seen that in the hniit A = one recovers the resuh obtained in [18]. 

One can also obtain the critical exponent relating A and Ac- For a specific 
value of Ac, one obtains two values for A. Denoting the difference of these values 
by A A, it is seen that 

AA^[Ae(0)-Ac]'/'. (32) 

3 The strong phase 

For A > Ac, there is a region where the density p^ exceeds one. Then pd cannot 
be pw and one uses another ansatz for pd, which is ps' 



Ps{z) 



1, z e [4,Cs] 

Ps(z), z £ [6s, ds] U [cs,as] 



(33) 



This is the strong phase. An argument similar to that used in [18] shows that 
A 



i/s(z) = ^{4z3 + (3z2)A 



^[{z - a,r - t2] [(z - s,)2 - i2] [4 (^^ + ss + z) + 3 A]} 



V[(^-^s)2-r2][(z-.,)2-t2] 



dC 



,^ {z-C)R{C)' 



(34) 



where 



s :- 



c + d 
c-d 



R{z) := ^[T^-{z-a,Y][tl-(z-s,Y] 



(35) 



Hs{z) should behave hke z ^ for z -^ oo. Expanding Hs{z) for large z, It is 
seen that 

Hs{z) ^aiz + ao + a^iz^'^ + 0{z^'^). (36) 

Hence one arrives at three equations 



One also has 



which is 



ai = 0, 

ao = 0, 

a_i ~ 1. 

dz[g{z)-H,{z)]^{), 



d. 



-\ dz[4(CTs + Ss + 2) + 3A]i?(z) + P 

d. 



Az 



rfe 



dC 



R{z) 



1 



,^ i?(C) z - c 



(37) 
(38) 

0. (39) 



Using this equation and the equations (37), one can obtain the four unknowns 
fs, Tg, Ss: and tg. Expanding these in terms of (A — Ac), and using (12), one 
can obtain an expansion for the derivative of the free energy with respect to the 
area. 

However, to obtain the behavior of the free energy in the strong phase at the 
limit A ^ 0, there exists a simpler way based on [18] and [27]. It is known from 
[27] that the transition is of third order, meaning that the difference between 
Hjyz) and H^{z) for large z behaves like {A — Ac)^: 



lim {[A-A,iX)]-^H,{z)-H„{z)]}^B 

A^A+{\) 



z, hm p^ 

A^A,{\) 



(40) 



Increasing A further, ps itself would exceeds one in sonic region. That is, a sec- 
ond phase transition occurs. Denoting the classical density after this transition 
by Ps 2 1 and the area corresponding to this transition by Ac 2 , one has 



lim {[A - Ac2(A)]-2 [H,^{z) - H^iz)]} = B 

A^A+^{X) 



z, lim Ps 

A^A,2W 



(41) 



Now consider the limit A ^ 0. At this limit the two transition occur at the 
same area Ac{0). One also has 



lim 



lim Ps 

A^Ao2(A) 



~ lim 



lim p^ 

A^A,(\) 



(42) 



So, adding (40) and (41) at the limit A ^ 0, one arrives at 



lim {[A-A,{0)]^ lim[Hs2{z)-H^{z)]} = 2B 

A— yl+(0) ^^0 



z, lim lim p^ 



Defining 



lim {[A AM? lim ^^(A) - F;(A)]}, 



it is then seen from (40) and (43) that 



lim {[A- AM? \immA)-FliA)]} = -I, 



A^AJiO) 



A^O 



(43) 
(44) 

(45) 



where (12) has been used the relate F' to H. 

Now consider the case A = 0, which has been investigated in [18]. When 
A = 0, the density p^ has two equal maxima and at the transition both maxima 
exceed one. It is then seen that [limA^o -P's 2(^)1 i^ ^ fact the derivative of 
the free energy in the strong phase for the case A = studied in [18]. So the 
difference between Fl{A) and i^4(^) m the hmit A ^- 0, is half what obtained 
in [18]: 



lim[^.'(A)-^;(A)] = ^ 



A- A, 



o 



A- A, 



(46) 



4 The second transition 

As it was pointed out in the previous section, p„ has two maxima which are 
equal when A = 0. If A is not zero, then one of these becomes larger and the 
absolute maximum of pw. The transition point is where this maximum is going 
to exceed one. After this transition, the classical density is ps, for which there 
is a region where the density is equal to one. However, in the region around 
the point the other maximum of p^ occurs, the density ps continues to increase 
and there may be some point where this second maximum exceeds one as well. 
Here the second transition occurs. For small values of A it is easy to find this 
point. Once again the result of [27] is used. There it is proved that the difference 
between the Ps{z) and Pw{z) vanishes faster than (A—Ac), provided the distance 
of z from the region where ps is one is large compared to the width of the region. 
This criterion is satisfied for z around the point where the second maximum of 
Pw occurs. So, 

Ps{z) = p^{z) + o[{A ^ Ac)], zKz^, (47) 

where Zm is the point where the second maximum of p„ occurs. This shows that 
up to leading order, the second transition occurs where the second maximum of 
Pw approaches one. So, 



Ar 



AM 



1 



AM 
47r 



|A|= 



(48) 



The phase picture of the system for small A is now complete. There are three 
phases: weak, strong, and stronger. 



• The weak phase. Here one has 



A < A,(0) 



Ac(O) 
An 



(49) 



the classical density is pw, and the classical density is everywhere less than 
one. 



• The strong phase. Here one has 



AM 



1- 



^c(O) |^|3' 



47r 



<A< Ae(0) 



1 



47r 



(50) 



the classical density is ps, and there is one interval where the classical 
density is equal to one. For fixed A, the area interval corresponding to 
this phase is seen to be proportional to A'^. This phase consists of two 
parts. For positive (negative) A the region where p^i is one is around z_ 



The stronger phase. Here one has 



^c(O) 



1 



^c(O) 3' 



<A, 



(51) 



the classical density is ps 2 , and there are two intervals where the classical 
density is equal to one. 

It is seen that if A vanishes, the strong phase disappears and only two phases 
remain, which agrees with the result of [18]. 

5 Concluding remarks 

A gYM2 with quartic and cubic couplings was studied. The effect of the cu- 
bic coupling on the classical density and the free energy was investigated. The 
quantitative changes in the weak phase were determined for small cubic cou- 
plings. The strong phases of the system were also investigated. It was seen that 
for small but nonvanishing cubic couplings there are two transitions, compared 
to a single transition if there is no cubic coupling. Both of these transitions are 
third order (in area). It was also seen that for small but nonvanishing cubic 
couplings, the jump in the third derivative of the free energy in each of these 
transitions is half the jump corresponding to the single transition which occurs 
if there is no cubic coupling. 
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